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Abstract 

We present a general study of a three neutrino flavour transition model based on the supersym¬ 
metric interactions which violate R-parity. These interactions induce flavour violating scattering 
reactions between solar matter and neutrinos. The model does not contain any vacuum mass or 
mixing angle for the first generation neutrino. Instead, the effective mixing in the first generation 
is induced via the new interactions. The model provides a natural interpretation of the atmospheric 
neutrino anomaly, and is consistent with reactor experiments. We determine all R-parity violating 
couplings which can contribute to the effective neutrino oscillations, and summarize the present lab¬ 
oratory bounds. Independent of the specific nature of the (supersymmetric) flavour violating model, 
the experimental data on the solar neutrino rates and the recoil electron energy spectrum are in¬ 
consistent with the theoretical predictions. The confidence level of the x^-analysis ranges between 
~ 10“"^ and ~ 10“®. The incompatibility, is due to the new SNO results, and excludes the present 
model. We conclude that a non-vanishing vacuum mixing angle for the first generation neutrino is 
necessary in our model. We expect this also to apply to the solutions based on other flavour violating 
interactions having constraints of the same order of magnitude. 


BAGS numbers: 12.60. Jv, 13.15.+g, i 4 . 6 O.Pq, i 4 . 8 O.Ly, 26.65.+t 


1 Introduction 


The most elegant solution of the solar neutrino problem []^ is the matter-enhanced neu¬ 
trino oscillation |^-[^ (for review articles see Hi)- It is based on the Mikheyev-Smirnov- 
Wolfenstein (MSW) resonance mechanism [Q which requires neutrino mixing in vac¬ 
uum. 

An attractive alternative interpretation of the solar neutrino dehcit is the neutrino 
flavour transition, due to the existence of some flavour changing (FC) and non-universal 
flavour diagonal (NFD) interactions between neutrinos and solar matter |^-[|^. In such 
a scenario, neutrino oscillations are induced by solar matter alone. Thus the neutrinos 
can even be massless and no neutrino mixing in vacuum is needed. 


* e-mail: dreinerOth. phy sik. uni-bonn. de 
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The explanation of the solar nentrino pnzzle in terms of FC and NFD scattering of 
nentrinos on solar matter, has been explored qnantitatively in |^-[^, and in for the 
particnlar case of nentrinos scattering on electrons characterized by the absence of a reso¬ 
nance. In all these stndies, the predicted nentrino rates are consistent with the resnlts from 
the solar nentrino experiments (Homestake, GALLEX, GNO, SAGE, SNO, Kamiokande 
and Snper-Kamiokande). Acceptable hts have also been fonnd in the framework where 
both matter-induced neutrino flavour transitions and neutrino mixing in vacuum exist 


12, 17 


There are also three flavour transition models, solving both the solar and atmospheric 
neutrino anomalies and in which FG and NFD interactions play a fundamental role. In 
TSl, 113 a model is constructed in which FG and NFD interactions exist between all 


neutrino flavours, all the vacuum mixing angles vanish and only the third neutrino is 
massive (in vacuum). The authors of have considered a scenario where FG and NFD 


interactions are present between the hrst and third neutrino flavours, only the vacuum 
mixing angle of — Vj. sector has a signihcant value and the second and third generation 
neutrinos are massive (in vacuum). These scenarios have been chosen so that there is a 
vacuum contribution to the mixing in the — Vr sector because the atmospheric neutrino 
problem cannot be entirely solved by FG and NFD interactions. The laboratory bounds 
on the interactions are too severe [^. The other main constraint, which these models 
satished, is the limit on the z/g — mixing angle from reactor experiments |^, 1^, ^ . 


In the present paper, we study the three flavour transition model where FG and NFD 
interactions are present between all neutrino flavours. However^ we restrict the vacuum 
neutrino masses and a mixing angle to the second and third generation. We thus address 
the question: do we need a vacuum neutrino mixing angle between the hrst and second 
generation to explain the solar neutrino problem or do FG and NFD interactions suffice? 
The model we study is the same model as in but with the most general FG and NFD 
interactions. 

We concentrate on the situation where the FG and NFD interactions are due to super- 
symmetry with broken R-parity [^, ^ , 


although our analysis is sufficiently general. 


that it can be easily applied to any new FG and NFD interaction which is of the same 
(small) order of magnitude. The superpotential is written in terms of the left-handed 
superhelds for the leptons (L), quarks {Q) and Higgs of hypercharge 1/2 (H) and the 
right-handed superhelds for the charged leptons (E'^), up and down type quarks (G'^, V^), 


Wm = 
R: 


= E 


+ K^,L,Q,DI + 


( 1 ) 


i, j, k are havour indices, Xijk, X'ijk, X'-jf. are dimensionless coupling constants and /ij dimen¬ 
sion one parameters. We consider all the relevant flp interactions, namely all the diherent 
havour conhgurations of both A and A' couplings, in contrast with the preliminary study 
p0| where only a small set of the A couplings were considered. 

Based on the most recent Standard Solar Model (SSM): “BP2000” p8[ (and including 


made in 


the new measurement of 517 ( 0 ) ^), we determine whether our generic model remains 

a possible interpretation to the solar neutrino anomaly. We have included the results 
from the Sudbury Neutrino Observatory (SNO) |^, |^, ^ In this sense, our work 
constitutes an update of the preliminary analysis in 
results from SNO appeared. 


20[| which was performed before the 
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2 Description of the Model 


2.1 Hamiltonian 


Our model is characterized by the time evolution equation of neutrino flavour eigenstates, 


.d 

i— 

dt 


I ^eit) 
V ^rit) 


H 


I ^e{t) \ 
V Mt) ) 


( 2 ) 


where the Hamiltonian has a vacuum contribution as well as a matter-induced part: 


H = 


E X l3x3 


( 0 
0 
0 




u 

^ cos 20-3 


V 


rn - - rn out) 

2 2 2 2 

mi—mi . r,/iv mi+mi , mi—mi r,nv 

Sin 20^3 cos 20^3 y 

^ -Rll -|- Hi -|- A 2 Ri2 Ri3 

Ri 2 R22 + H2 R23 

\ -R13 R23 R33 + H2 


(3) 


E is the neutrino energy, I 3 X 3 is the 3x3 identity matrix, and m| are, respectively, 
the second and third generation neutrino masses (in vacuum). 0^3 is the vacuum mixing 
angle in the — Uj. sector. Hi and H 2 are the contributions due respectively to the 
and boson exchanges. Note, that we have omitted any hrst generation contributions 
to /{2 E))yacuum, as discussed in the introduction. 

The Rij terms arise from FC (if Rij ^ 0 [i ^ j]) and NFD (if i?ii 7 ^ R 22 7 ^ R 33 ) 
interactions predicted by new physics: here R-parity violating supersymmetry. Hi = 
y/^GpUe and H 2 = —GpriN/V^ and Gp is the Fermi coupling constant, ng is the electron 
density and np; the neutron density in the sun. Since we concentrate on ftp interactions, 
the Rij terms are given by 


Rij 


Rij{X) 


^p(A') 


Rij{X) + Rij{X') , 
XiklXjkl 

4m2([ 


-n. 


E 

I 


XiiiX 


ju 




n. 


E 


/ A' A' 
\4:m?{dm) 


Hd H- ^ - Ud 

4m2(d^) 


(4) 


where Ud is the down quark density and {m?{dm)) is the squared mass of the 

slepton If I (squark dm) exchanged in the scattering reaction Ui+e {r'i+d Uj+d). 

The scattering reaction also hxes at least one of the generation indices to be 1. The second 
term of Rij{X) vanishes if i = 1 or j = 1 , due to the anti-symmetry of Xijk in the hrst 
two indices. In the following, we perform a general analysis with the Rij. Given the 
approximations discussed in Section pH| , this analysis can be applied to any new FC and 
NFD interactions. 
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2.2 Diagonalization 


Since the diagonal contribntions to the Hamiltonian do not affect the ffavonr transition 
mechanism, only the following effective mass sqnared matrix, 




— H — (^E + Ai + A2) X 1 


3x3? 


is relevant. It can be diagonalized as. 


(5) 


( 6 ) 


The unitary matrix U is parameterized by the complex phase 6 and the mixing angles 
6 ^ 12 , 9 i 3 and 623 as. 



1 

( 

~ 2 

0 

0 \ 

— 

H7= — 


0 

rhn 

0 

\2 Ea 

2E 

\ 

0 

Z 

0 

m\ J 



( 1 

0 

0 \ 

/ 

Cl3 

0 

si3e \ 

/ 

C12 

S12 


U = V 23 X V13 X V12 = 

0 

C23 

S23 


0 

1 

0 


-S12 

C12 



lo 

—"S23 

C23 / 

1 

-8x3^^ 

0 

Cl3 / 

1 

0 

0 

1 / 


/ 


U = 


C 12 C 13 


— S 12 C 23 — Ci2S23'Sl3e 
ifi 

\ S12S23 - Ci 2 C 23 Si 3 e“ 


,i5 


’S 12 C 13 

C 12 C 23 — Si2S23Sl3e*' 
— C12S23 — Si 2 C 23 Si 3 e 


iS 


sise 

'S 23 C 13 

C 23 C 13 


(7) 


Here Sij = sin 6 ij and = cos 6 ij. In the following, we choose the parameters of the 
Hamiltonian H to be real, so that 5 = 0 above. First, we rotate the effective mass 
squared matrix by V 23 x V 13 and obtain in the new basis. 


/M^\ 

f Rnch 

— a sin 2013 + II+S 13 

/3ci3 

0 \ 



f3ci3 


/3si3 

\ 

0 

Psi3 

i?iiSi 3 + a sin 20 i 3 + H+C 13 / 


with. 


H 4 - = 


a = i?l2-S23 + 7 ?i3C23, (3 = R 12 C 23 — 7?13S23, 

_ + ^2 A I 7?33 + R 22 

= i -7T-Hi + - 


( 8 ) 


(9) 


2 cos 2(0,^3-023) 

^- - -cos 2023 + R 23 sin 2023 ]. 


( 10 ) 


The mixing angles of Eq. i 

tan 2023 = 

where Am^, = — m,, and. 


(p)0|) are given in terms of the Hamiltonian parameters 
sin 29 23 ^^. 77112 /‘2'E + 2 R 23 


cos 20,3 Am 32 / 2 £' + R 33 — R 22 


( 11 ) 


tan 2013 = 


2 a 


n.-R 


( 12 ) 
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We show at the end of Section |2.3| that, in our framework, the mixing angle of the 
Ue — Vr sector has a negligible value 


0 


13 


0 


(13) 
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Therefore, the upper left 2 x 2 part of matrix (j^) can be readily diagonalized via the 
mixing angle 6*12 dehned by, 


tan 26*12 ~ 


_2^Ci3_ 

— -R11C13 + a sin 26*13 — f^+^is 


Then, the effective mass eigenvalues are. 


( 14 ) 


Ri (/?iiCi 3 -asin 26 * 13 + fl+Si 3 )ci 2 -fesin 26*12 + fi_Si2 , ( 15 ) 

Ztj 

~ 2 

Tin 

^ ^ (/?iiCi3-a sin 26 * 13 + fl+Si3)si2 + fe sin 26*12 + fl-Ci2 , ( 16 ) 

Ztj 

~ 2 

TTT 

^ ^ i?iiSi3 + asin26*i3 + fl+Ci3 . (17) 

Ztj 

At this stage, an important comment must be made concerning the rnixing 

angle 6*12. We see from Eqs. (|]) and (0) that tan 26*12 is proportional to the sum of two 
terms -R12C23 and —-R13S23. R12 constitutes the off diagonal contribution to the Ve — 
sector of the Hamiltonian in Eq. (j^). The second term results from a transmission (due to 
i?i3) of the mixing in the sector (dehned by the angle 623) into the — z/^ sector. 

We conclude that the i?i2 and R13 contributions, induced by interactions, generate an 
effective mixing in the Uf, — z/^ sector and thus play a fundamental role in the neutrino 
havour transition mechanism. 


2.3 Approximations and Constraints from Reactor and Atmospheric Neu¬ 
trino Experiments 

The less stringent bounds on the products of ftp coupling constants are typically of order 


35-39 


A-A ~ 10 “^ 


m 


f 


100 GeV 


A 6 {A, A', A"} , 


( 18 ) 


where represents the typical scalar superpartner mass. Hence, Rij is at most of order 
(c./. Section for a more precise discussion): 


Rij ^ Hx X 2.5 ■ 10 ^GeV , na,e{ne,nd}. 


rid, fhe maximal electron density reached is 


(19) 
4.6 ■ 


Moreover, inside the sun, n, 

10 ^^ eV^ 12 ^ and the energy of the produced neutrinos is never higher than E = 20 MeV. 
Therefore, 


Ri 


Am 32 / 2 E 


< (2(10 


-S', 


( 20 ) 


since Am'^2 ~ 10 ^ eV^, as we will see in Eq. (P 7 |). In the following, we perform an analysis 
of the Hamiltonian (^ with general entries Rij, however with the restriction ( 0 ). above. 
Thus our analysis applies to any new physics contributing to the Rij, which respects the 
bound 0 ), t.e. with couplings obeying the bounds (B- 

We deduce from Eq. (O) that the experimental data on atmospheric neutrinos M, H 


are consistent with the Up —>■ Vj. oscillation scenario based on an effective two flavour mass 
matrix, for the following values. 


Am32 = rug — 11X2 ~ 10 eV 


( 21 ) 
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and, 


TT 


^23 


( 22 ) 


Eqs. (pd]), pop and reqnire for the vacnnm mixing angle 6^23, 


0, 


23 


TT 

~ ^ ~ 023 


( 23 ) 


We also see from Eq. (|T^ that the ratio of Rij and Ai reaches at most the value, 


( 24 ) 


From Eqs. (|^, (p^), 
approximation as. 


and 


^ < 10-2 . 

4 li 

it is clear that the mixing angle 0i2 is given to a good 


tan 2012 


\/2{Ri2 — i^is) 


2E 




( 25 ) 


From Eq. (PB| ) we see there is potentially a resonant point at 0i2 = vr/d. We hnd that 
the associated resonance condition, namely m^l2E = Ai, can be indeed satished inside 
the snn. However, for the atmospheric nentrinos, which have energies of Cl(GeV), this 
resonance condition cannot be fnlhlled. The relevant parameter valnes are: m\ ~ 5 ■ 
10 “®eV^ (see Section |^) and the electron density inside the earth nf ~ (3 — 6) cm“^ A/^, 
where Ma is Avogadro's constant. 

From Eqs. (H), dD, (0), m and (p^), we derive the following approximation for 
the effective mass sqnared difference Aml2 (c./. Eq. (pT|) ), 


Am 


32 


Am 32 


We then dednce from Eqs. (pp and (| 26 |) that, 

Am32 ~ 10 “^ eV 

Besides, Eqs. 


( 26 ) 


( 27 ) 


, (HD, (PD, 


and 


Am2i = ml 


~ 2 2 
m-i ~ nrir, 


imply that, 

2EAi. 

Finally, we jnstify the approximation of Eq. 
in Section First, from Eqs. (p 3 |) , (p^ and 
013 : 

y/2{Ri2 + i?i3 


( 28 ) 

, nsed both in the present section and 
, we obtain a good approximation for 


tan 2013 


Ai _ 4, 


( 29 ) 


2E 

Now, we see from this resnlt and Eq. (p^) that 0i3 ~ 0 ontside of the resonance. The res¬ 
onance point associated with 0i3 cannot be reached inside the snn, becanse the resonance 
condition m\/ 2 E = Ai cannot be satished. Indeed, inside the snn (c./. Eq. 


mt 


> 


Am| 


32 > 


2.5 ■ 10 "^^ eV , and, Hi ~ 7.5 10 "^^ eV 


-12. 


( 30 ) 


2 E 2 E 

The resonance condition m‘l/ 2 E = Hi can also not be fnlhlled for atmospheric nentrinos, 
given nf above. Moreover, it is clear from Eq. (P^ that in vacnnm, the 0i3 mixing 
angle is exactly eqnal to zero (like 0i2), consistent with the constraint obtained at reactor 
experiments (with an ehective two havonr mass matrix): sin^ 20i3 ~ 0.1 for ^m\i ~ 
3 ■ 10 “^eV^ [p 2 , M]. Therefore, within the present framework, 0i3 approximately 

vanishes from the point of view of the reactor, solar and atmospheric nentrino experiments. 
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3 Solar Neutrino Rates 


3.1 Theoretical Transition Probabilities 


transitions of solar neutrinos 


Since to a good approximation 6*13 is vanishing, the Ue 
obey a two flavour transition dynamics described by the upper left 2 x 2 block of the 
matrix in Eq. (H). Thus, in the non-adiabatic transition case, the solar z/g —> transition 

probabilities = 1 — are dehned by |^, ^ 


Pv.^u, - ^ (1 - (1 - 2Px) cos^ 20*2 I , 


Px > 


(31) 


for a neutrino produced above the resonance density, and — 1 , for a neutrino 

produced below the resonance density. In the adiabatic transition case [0], 


Pue^ue - ^ ( 1 + (1 - 2Px) cos 20^2 cos 20 


Px<- 

X 2 


(32) 


In Eq. (^p and Eq. (^), 0 i 2 and 0*2 represent the matter-induced 0 i 2 mixing angle at 
the neutrino production point and solar surface, respectively. P^ is the Landau-Zener- 


Stueckelberg transition probability [44 


Pr ~ 


(33) 


In this equation, P ~ 1 — tan^ Q, and 0]’“ is the 0 i 2 mixing angle in vacuum, since 
the electron density inside the sun is approximately given by an exponential function of 
the distance r to the solar center [^. From Eq. ( |2^) we deduce that 0^1'^ ~ 0, so that 
F ~ 1 here. The 7 parameter of Eq. (P3|) is given at the resonance point. 


7 = 


Am2i/2E 


o d0i2 
dr 


(34) 


The resonance point is determined by (c./. Eq. (|^)), 


— = 44i. 

2E 


(35) 


Let us derive the expression for 7 within our model. At the resonance point, 012 = 7r/4 
so that Eqs. (P^ and (P^), together with Eqs. (0 and (P3|), give. 


Am|i 


2E 


= V2\R 


12 


R 


13 Ires 


(36) 


From Eq. ([25| ) we see that |2 dOu/drlres = GF\^/iRi 2 — Pi 3 )|res and thus. 


7 = 


V2\R 


12 


R 


131 


Gf 


drie I 
dr 


(37) 


3.2 Fit with Experimental Results 

In this subsection, we present the free parameter values giving rise to the highest degree 
of consistency between our theoretical predictions and the experimental solar neutrino 
rates. The relevant free parameters in Eqs. (^,|^) are and R, where 


(38) 


|fli2(A)-fii3(A)| 


rig 
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for A-/?p interactions and 


R = 


|i?i2(A')-i?13(A')| 


nd 


(39) 


for X'-^p interactions (c./. Eq. (|[)). In certain specified cases, we also consider fs as a 
free parameter: fs is defined by = /b x ^ where is the flux of solar electron- 
neutrinos produced in the decay ^ ^Be* + e~^ + u, 
prediction within the SSM [ 


and ^8^^ is the theoretical 


30|. This is equivalent to allowing for an arbitrary 


normalization of the <l>8s flux. 

In Table |^, we show the values of the free parameters for which the best fit is obtained 
between the theoretical solar neutrino rates defined as, 

dgi{E) dhi{r) 


= dE dr 


dE dr 




(40) 


and the corresponding present experimental solar neutrino rates P5||-|^ 


The index j runs over the 5 solar neutrino experiments (we use the combined result for 
the GALLEX, GNO and SAGE experiments |^): j ={Ilomestake, GALLEX -|- GNO -|- 
SAGE, SNO, Kamiokande, Super-Kamiokande}. 

Let us briefly describe the parameters entering Eq. (^). First, is the theoretical 
prediction for the solar neutrino detection rate in experiment j, if neutrinos have a flavour 
transition average probability Piy^^y^{E, r). The (pi are the flux of solar electron-neutrinos, 
where i labels the source reaction: i ={pp,^Be, pep, ^^0, ®B, hep}. dgi{E)/dE 

and dhi{r)/dr respectively represent the neutrino energy and the initial radial position 
distributions (normalized to unity). Both these fluxes and distributions are taken from 
the most recent SSM (BP2000) @-|^. Finally, the (Tj{E) denote the neutrino detection 
cross sections corresponding to the different solar neutrino experiments |^, . 

In the case of the Kamiokande and Super-Kamiokande experiments, the product crj{E) 
Py^^t,XE,r) in Eq. (EO) stands for. 


'Tth 


dT / dT' e(T') p(T, T') 


d(Tf{E,T') 


da]'^{E,T') 

dT' 


dT 

l-K 


P. 


{E,r) 


(41) 


where T' is the true energy of the recoil electron, T its value measured experimentally, 
Tth the energy threshold of the experiment, T^ax the maximal energy of the experiment, 
e(r') the detection efficiency, p(T,T') the energy resolution function and da'^''{E,T')/dT' 
{daj'^{E,T')/dT') the differential cross section for the solar electron-neutrino (muon- 
neutrino) elastic scattering reaction on electrons |l[]. The resolution function of Eq. (^I[) 


can be satisfactorily represented by the following Gaussian function 


p(r,r') = 


-{t-t'P 

g 2A^(T') 


The standard deviation of the energy measurement is approximately given by. 


A (TO 


/lOMeKy/2 


(42) 


(43) 


We take Tth = 7.5 MeV T^ax = 20 MeV 


for the Kamiokande detector, and Tth = 5 MeV |50 


i 0.7 0 and a = 0.143 ||, |5|, |5|, |5l, 
Tmax = 20 MeV 13, e 0.7 0, |58 


and a = 0.153 |59, 60|, for the Super-Kamiokande one. 
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For the SNO experiment, the product aj{E)Py^^^^{E,r) in Eq. (|i0|) stands for, 

'-pSNO 

I m. n. T. 


frpSNO 

-^th 


/■oo ^daf^{E,r)- 

dT dT' p{T, T') ^ 


e — 


(E,r). 


( 44 ) 


daf^{E,r)/dr is the differential cross section of the charged current reaction between 


the solar electron-neutrinos and deuterium: i>ePD^p + p + e [53]. The energy 


resolution function p{T,T') of Eq. (^) is dehned as before by Eq. (|4^), but with the 
following standard deviation A(T') 


A(T') = - 0.4620 + 0.5470 ^T'/lMeV + 0.008722 (T'/l MeV) 


MeV. 


(45) 


In the present analysis, we use = 6.75 MeV -|- rrie 

= 13 MeV + me 


, me is the electron mass, and 


31 


The points of parameter space corresponding to the best £t between the theoretical 
and experimental solar neutrino rates have been obtained by performing a y^-analysis. 
Following the analyses of |^, we dehne our y^-function as. 


xi=E (d'* - fi:'‘)\o%pps - 'P.D. 


(46) 


n,j2 


where ji and j 2 run over the solar neutrino experiments and is the squared error 

matrix containing both the experimental (systematic and statistical) [^I], ^ 


50| and theoretical errors (computed according to the analyses of Q) on the solar 
neutrino rates. The theoretical errors depend on the uncertainties of both the neutrino 


fluxes P 8 


68 [| and the detection cross sections [^ 


In Table ^ our results concerning the fits of solar neutrino rates are given in two 
different contexts: before (case ( 1 )) and after the publication of experimental data from 
SNO (case (2)). In the “post-SNO” context, we include into our yVanalysis the 


rate of solar electron-neutrinos detected at the SNO experiment via the charged current 
reaction on deuterium |^. Furthermore, within this context, we £x the normalization 
factor Jb to unity, following the study of |]^. The reason is that the value of the neutrino 
flux $85 measured at the SNO experiment (with the flux of neutrinos detected through 
the neutral current reaction on deuterium) [^, |^, is in agreement with its prediction 
in the SSM ||2^, |^: Finally, let us note that within this post-SNO context, 

we do not implement into the y^-analysis the rate of solar neutrinos detected at the SNO 


experiment via the elastic scattering reaction on electrons ||31| , M . This is justihed by 
the fact that this rate is consistent with the rates of neutrinos detected at the Kamiokande 
and Super-Kamiokande experiments via the same scattering reaction Q, which are 
already taken into account in our analysis. 

In Table |I|, we also present the values of the parameters corresponding to the best £t of 
all 5 solar neutrino rates and the recoil electron energy spectrum (case (3)). This £t takes 
into account the consistency between the recoil electron energy spectrum obtained with 
the Super-Kamiokande detector from 1258 days of data [^ and its theoretical prediction. 
These results have also been derived through the y^-method with the y^-function 


xl 






(47) 


31,J2 


Here the [j = 1,..., 19] are the solar neutrino rates observed in the Super-Kamiokande 


experiment corresponding to the 19 bins of a given recoil electron energy range [50] and 
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Model 

(fit) 

mi 

(10-5 

eW) 

R 

(10-24 

eE-2) 

fBla 

r 

{DOF; C.L.) 

H 

(0.297 

±0.026) 

G±G±S 

(0.557 

±0.039) 

SNO 

(0.295 

±0.023) 

K 

(0.472 

±0.064) 

SK 

(0.391 

±0.014) 

A 

(1) 

0.545 

1.319 

Jb = 
1.087 

0.19 

(1; 66.3%) 

0.295 

0.574 

X 

0.476 

0.390 

A 

(2) 

0.625 

1.055 


6.25 

(3; 10.0%) 

0.244 

0.595 

0.316 

0.402 

0.369 

A 

(3) 

0.648 

0.561 

a = 

0.624 

46.60 

(21; 0.11%) 

0.560 

0.645 

0.704 

0.738 

0.713 

A' 

(1) 

0.534 

1.112 

/b = 
1.115 

0.18 

(1; 67.1%) 

0.303 

0.562 

X 

0.491 

0.389 

A' 

(2) 

0.625 

0.814 


6.23 

(3; 10.1%) 

0.244 

0.595 

0.316 

0.402 

0.369 

A' 

(3) 

0.386 

0.257 

a = 

0.459 

51.73 

(21; 0.021%) 

0.865 

0.629 

0.997 

0.992 

0.972 


Table 1: The values of the parameters m^, R, Jb and a (see text) associated with the best fits between 
the solar neutrino experimental data and the theoretical prediction. They are given in the 3 cases we 
consider: (1) the solar neutrino rates except those of the SNO experiment, (2) all the solar neutrino rates, 
including the rates detected at the SNO experiment via the charged current reaction, and (3) all the solar 
neutrino rates together with the recoil electron energy spectrum obtained at the Super-Kamiokande 
experiment. Both the models based on the interactions of type A and X' are considered. We show 
the values, computed with the obtained best-fit parameters, of the relevant y^-function (with associated 
number of degrees of freedom denoted as DOF and confidence level denoted as C.L.), namely either Xb 
(for cases 1 and 2) or y| (for case 3). We also indicate, for the best-fit parameters, the theoretical solar 
neutrino rates (see Eq. (|^) normalized to their value expected in the SSM BP2000 at 

the considered experiments: Homestake (H), GALLEX-I-GNO-I-SAGE (G-l-G-l-S), SNO, Kamiokande (K) 
and Super-Kamiokande (SK). The corresponding present experimental solar neutrino rates (together 
with their experimental uncertainty) are written in parentheses. Finally, we note that 

when the value of Jb or a (which are presented in the same column) is not given, it means that it is fixed 
to unity and/or it is not relevant. 


the [i = 1,..., 19] are their predicted values. The are calculated with Eqs. 

I) (applied to the Super-Kamiokande experiment) by integrating T over the relevant 19 
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energy intervals 


The matrix \a 


Sijlj2 


entering Eq. (^71) is equal to |T4 


^jlj2iWstat\jlj2 T Wuncorr\jlj2) T [(^corr]jlj2 T [^th\jlj2 


(48) 


The [a 


stat\ 3 i 32 


represent the squared experimental statistical errors 0, 


Wcorr\hj 2 and 




uncorr\3i32 


are respectively the correlated and uncorrelated contributions to the squared 

Wth]jij 2 are the squared theoretical uncertainties 


experimental systematic errors |0 and 


|H^ . Finally, the overall normalization factor a introduced in Eq. ( ^71) has 
been taken as an additional free parameter in the £t, in order to avoid double-counting 
with the data from the Super-Kamiokande experiment on the total event rate which is 
already included in Xr- 

Let us comment on the results of the y^-analyses given in Table |I|. In case one considers 
the solar neutrino rates excluding the SNO experiment, case (1), one obtains acceptable 
hts between the theoretical and the experimental rates, namely hts at 66.3% and 67.1% 
of Conhdence Level (C.L.) in the scenarios containing A and A' interactions, respectively. 
When adding the SNO results to the analysis, the £t grows signihcantly worse: the theo¬ 
retical predictions are now compatible with the experimental data only at the 10.0% C.L. 
and 10.1% C.L., respectively. Finally, considering all the solar neutrino rates (including 
the SNO data) as well as the recoil electron energy spectrum, case (3), the obtained hts 
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have a confidence level of 1.1 x 10 ^ C.L. and 2.1 X 10 C.L., respectively: the theo¬ 

retical and experimental results are not compatible in this case. This incompatibility is 
due to the SNO results and not the recoil electron energy spectrum. In the authors 
obtained an acceptable fit, without the SNO data but including the recoil electron energy 
spectrum. 

Furthermore, we observe in Table Q that the values of the quantity R, corresponding 
to the best fit, are smaller in case the solar neutrino flavour transitions are induced by 
ftp interactions of the type A' instead of A. The reason is that the quantity {c.f. Eq. 

= I-R12 ~ -Rial = |-Ri 2(A) — i?i3(A)| = RUf, , 


(49) 


entering the transition probability (see Section |3.1|) , in case neutrino flavour tran¬ 

sitions are due to A couplings, is replaced by {c.f. Eq. (|3^)), 

= |Ri 2 — -Rial = |-Ri 2(A0 — -Ria(-^0I = -R^d • 


The decisive difference is that 0.50 ~ Ue ~ 0.78 inside the sun ^ 

words, a modification of the matter-induced transition probability Pi,^- 


(50) 

In other 
due to an in¬ 
crease of the solar matter density can be compensated by a decrease of effective neutrino 
ftp interaction strength (quantified here by R). Let us note that this compensation is not 
exact, as the solar matter density depends on the radial position, while this is not true for 
the parameter R. This is why, in Table |^, the best-fit values of the other free parameters, 
and fs (or a), are also different in the two situations where flavour transitions are 
due to A and X' couplings. 

Let us make a brief final remark: we have found that, in the scenario where solar 
neutrino flavour transitions are simultaneously induced by ^p interactions of type A and 
A', the fits between observed and predicted solar neutrino rates and recoil electron energy 
spectrum are not significantly improved compared to the fits for which the results are 
shown in Table |l|. 


3.3 Comparison of the Results with the Bounds on the flp 
stants 


Coupling Con- 


In this section, we check that the relevant present limits on flp coupling constants are 
compatible with the values of R giving rise to the best fits of solar neutrino rates (see 
previous section). This check is also a self-consistency check. It allows us to verify that the 
typical bound on the relevant ftp coupling constants ( 0 ), on which the approximations 
(see Section |2.3|) used in the mass matrix diagonalization are based, are correct. 

• A coupling constants: For A-couplings R reads explicitly {c.f. Eq. (^)), 


R = 


A131A231 


Ai91 a 


121/^321 


4m2(f=*=) 4m2(/i^) 


(51) 


The largest best fit value for R in Table I is 1.319 10 The bounds on the R-parity 

violating couplings are usually given as bounds on a single coupling or on a product of 
two couplings. This is based on the assumption that the other couplings vanishQ. We 
adopt this approximation in the following. 

No bound exists on the ftp coupling constant product A 121 A 321 . The present exper¬ 
imental constraints at 2 cr on the single fip coupling constants A 121 and A 321 are given 

^The assumption that a single or a pair of p.p coupling constants are dominant is often adopted in the literature, as 
a simplification. This hypothesis is justified by the analogy between structures of and Higgs Yukawa interactions |70| . 
The latter exhibit a strong hierarchy in flavour space. 
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respectively by |^, |38|, p9|, [7I|] , 


Ai2i < 0.049 


m 


(ei) 


100 GeV 


A 321 < 0.07 ^ 


m[e 


R) 


100 GeV 


(62) 


Note that a different slepton mass appears in R and in the bound above. If we assume 
that m{jl) = we obtain from the second term in Eq. (pl|) that i? < 9 ■ 10 eV ^ 

and this model irrelevant for EG and NED interactions in the sun. If however m{fl) = 
m(eij)/4, which could very well be the case with non-universal boundary conditions at 
the unification scale, 

1.372 ■10■2^eV■^ (53) 

and the A 121 A 321 interactions can reach the best fit values for R. 

The strongest bound on the Yukawa coupling constant product A 131 A 231 is [[T^, 


A131A231 < 6.6 10 


-7 


m{UrL) Y 
lOOGeVy ■ 


(54) 


In this case we would need to require m{i'rL) ~ 300 x m{f) in order to obtain a sufficiently 
large R. For m(f) ~ 100 GeV (the LEP bound |^), this is not compatible with a 
supersymmetric solution to the hierarchy problem. We conclude that these interactions 
can not significantly contribute to EG and NED interactions in the sun. 

In summary, only when R is given by i? = |Ai2iA32i|/4m^(/l^), the present constraints 
on relevant fip interactions allow R to take the relevant value of 1.319 10“^"^eV“^. 

A' coupling constants: Let us now consider the model (not studied in paper p0| ) 


in which neutrino flavour transitions are due to A' ftp -interactions via solar down quarks. 


R = 


E 


/ \' \' V V V V V V \ 

I ^llm^21m _ ^lml^3ml _ ^llm^31m ) 


4md{d^)) 


(55) 


Based on the present constraints on fip interactions, we determine whether R can reach 
1.11 which is the maximal value appearing in Table 2. As before, we consider 

separately each term of Eq. (|55D, supposing that it dominates over all the others. 


The product A)^ 2 i'^ 32 i (third term in R with m = 2) is only constrained by the present 


individual experimental limits at 2 a on A) 2 i and A 321 


A;2i < 0.043 


m{dR) 

100 GeV 


A 321 < 0.52 


m{dR) 

100 GeV 


(56) 


Again a different squark mass appears in i?, m{s), than in the bounds, m{dR). If we 
assume they are equal, then 

5.6-IQ-^^eV^ (57) 

This is just below the maximal best-fit value. For a moderate non-degeneracy of m{dR) = 
1.5 X m{s), the maximal best-fit value can be obtained. 

The Yukawa coupling constant product Ajgj^Aggj^ (third term in R, with m = 3) is only 
constrained by the present individual experimental bounds on Ajgj^ (at 3 (t) [^, and 
on A^ 3 i (at 2 a) [||, [ 


74 


(10015^ 


A 331 < 0.45 at m{d) = 100 GeV . 


Assuming all squark masses are equal 

R<2.1- IQ-^^eV^ 


(58) 


(59) 
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which is below the best-£t values in Table |^. Given a non-degeneracy of m{d) = m{tL) = 
2.3 X m{b), the fit values for R can be reached. 

For m = 1 we have for the last two terms the coupling which has the strict bound 


at 2a pa, 75, 76 


A'ln < 5.2 10 


_4 f m{e) 


VlOOGeV; VlOOGeV 


mix' 


Oi \l/2 


(60) 


In these cases, we can therefore not reach the £t values for R, unless we require an extreme 
hierarchy in superpartner masses. 

The present experimental bounds at 2a on A 342 and A 442 (fourth term in R 

with m = 2 ) are. 


-^312 < 


o il ( Min) ^ 1 ,, ^ 0,021 ( ) , (61) 

UOOGeVj ’ \100GeVj ^ ’ 

In this case we always have the same squark mass and we can deduce a rigorous bound 

5.8-10“^®eV% (62) 


which excludes any relevant contribution to FG and NFD interactions in the sun. 
In the same way, the present experimental limits at 2a on A 343 and A 443 [^, p^, 


(63) 


result in i? < 5.8 ■ 10“^® eV“^ , precluding any signihcant contribution. 

The three products of Yukawa coupling constants ,m = 1,2,3, (second term 

in R) suffer the following strong constraint |78| , 


A' A' <17 in-7 f ^idkR) Y 

An,A2i.<l-7 10 


(64) 


Hence, R < 4.25 ■ 10 ^°eV ^ and no interesting solutions exist. 

The products ,m = 1, 2, 3 (hrst term in R) also have a stringent limit |]^ 


Kn^2n < 1-6 10 


-7 ( mjujL) y 

\W0GeV) 


(65) 


precluding any signihcant contributions to the FG and NFD interactions in the sun. 

In summary, if the two coupling constants A '^21 ^^0 -^3211 -^131 -^3311 dominant 

so that the quantity R is given by i? = |A32iA32i|/4m^(s), or R = |A^ 33 A 334 |/ 4 m^(&), 
respectively, the present bounds on relevant coupling constants allow R to reach the 
relevant value: 1.112 10 “^^eIA“^. 


4 Conclusions 

We have studied a three neutrino Favour transition model (see Section ^T|) , based on the 
supersymmetric /Zp—interactions between neutrinos and solar matter (see Section p.2D . 
We have not included a vacuum mass or any vacuum mixing angles for the hrst generation 
neutrino. Instead the ehective mixing is induced via the Favour violating interactions. 

Within this scenario, we have found that, among all the (/?p ) contributions Rij to the 
matter-induced part of the ehective mass matrix in Eq.(pD, only R 12 and i?i 3 play a role 
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(see Section |2.3| ). Then we have shown, from a systematic study of the constraints on 
the ftp couplings, that only the products A 121 A 321 , A( 2 iA 32 i and which enter i?i 3 , 

can result in a sufficiently large value of R in agreement with the best ht obtained for the 
solar neutrino data (see Section |3.3| ). 

However, whatever the flavour conhgurations of the ftp contributions are, we have 
found that the best hts of all solar neutrino rates together with the recoil electron energy 
spectrum have a conhdence level ranging between ~ 10“^ (scenario with A couplings) 
and ~ 10“^ (scenario with A' couplings). The discrepancy between the experimental re¬ 
sults and the theoretical predictions is due to the results from the SNO experiment. The 
considered model, which provides a realistic interpretation of the atmospheric neutrino 
anomaly and respects the constraints on neutrino oscillations obtained at reactor experi¬ 
ments (see Section |0|) , does not constitute an acceptable solution to the solar neutrino 
problem. This analysis does not rely on the nature of the /?p—interactions and thus 
applies to any new physics contributing to R 12 and i?i 3 , respecting the upper bound in 
Eq. ([T^). We conclude that in this case a non-vanishing vacuum mixing angle for the hrst 
generation neutrino is necessary. A future interesting investigation would be to determine 
whether the same model, but based on other types of lepton flavour violating interactions 
underlying the Standard Model, can represent a reasonable solution to the atmospheric 
and solar neutrino puzzle, or whether all such models can be excluded. 
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